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Abstract 

A probabilistic framework is introduced that represents stylized banking networks 
and aims to predict the size of contagion events. In contrast to previous work on 
random financial networks, which assumes independent connections between banks, 
the possibility of disassortative edge probabilities (an above average tendency for 
small banks to link to large banks) is explicitly incorporated. We give a probabilistic 
analysis of the default cascade triggered by shocking the network. Wc find that 
the cascade can be understood as an explicit iterated mapping on a set of edge 
probabilities that converges to a fixed point. A cascade condition is derived that 
characterizes whether or not an infinitesimal shock to the network can grow to a finite 
size cascade, in analogy to the basic reproduction number Rq in epidemic modeling. It 
provides an easily computed measure of the systemic risk inherent in a given banking 
network topology. An analytic formula is given for the frequency of global cascades, 
derived from percolation theory on the random network. Two simple examples are 
used to demonstrate that edge-assortativity can have a strong effect on the level of 
systemic risk as measured by the cascade condition. Although the analytical methods 
are derived for infinite networks, large-scale Monte Carlo simulations are presented 
that demonstrate the applicability of the results to finite-sized networks. Finally, we 
propose a simple graph theoretic quantity, which we call "graph-assortativity" , that 
seems to best capture systemic risk. 

Key words: Systemic risk, banking network, contagion, random graph, cascade 
condition, credit risk, financial mathematics, assortativity. 
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1 Introduction 



The study of contagion in financial systems is very topical in light of the recent global 
credit crisis and the resultant damage inflicted on financial institutions. "Contagion" 
refers to the spread of defaults through a system of financial institutions, with each 
successive default causing increasing pressure on the remaining components of the 
system. The term "systemic risk" refers to the contagion-induced threat to the fi- 
nancial system as a whole, due to the default of one (or more) of its component 
institutions. 

It is widely held that financial systems (see Upper 2011 and references therein), 
defined for example as the collection of banks and financial institutions in a developed 
country, can be modelled as a random network of nodes or vertices with stylized bal- 
ance sheets, connected by directed links or edges that represent exposures ("interbank 
loans"), each edge with a positive weight that represents the size of the exposure. If 
ever a node becomes "insolvent" and ceases to operate as a bank, it will create bal- 
ance sheet shocks to other nodes, creating the potential of chains of insolvency that 
we will call "default cascades". Financial networks are difficult to observe because 
interbank data is often not publicly available, but studies have indicated that they 
share characteristics of other types of technological and social networks, such as the 
world wide web and Facebook. For example, the degree distributions of financial 
networks are thought to be "fat-tailed" since a significant number of banks are very 
highly connected. A less studied feature observed in financial networks (and as it 



happens, also the world wide web) is that they are highly "disassortative" (see New- 



man 



[2002] ). This refers to the property that any bank's counterparties (i.e. their 
graph neighbours) have a tendency to be banks of an opposite character. For exam- 
ple, it is observed that small banks tend to link preferentially to large banks rather 
than other small banks. Commonly, social networks are observed to be assortative 
rather than disassortative. Structural characteristics such as degree distribution and 
assortativity are felt to be highly relevant to the propagation of contagion in networks 
but the nature of such relationships is far from clear. 

Our aim here is to develop a mathematical framework that will be able to de- 
termine the systemic susceptibility in a rich class of infinite random network models 
with enough fiexibility to include the most important structural characteristics of real 
financial networks, in particular with general degree distributions and a prescribed 
degree of edge-assortativity. Our starting point will be the Gai and Kapadia (here- 
after referred to as GK for short) cascade model of Gai and Kapadia [2010 and the 
analytical methods developed there and in Gleeson et al. 2011| for this model. The 
basic assumptions introduced in the GK model are: 

1. The network is a large (actually infinite) random directed graph with a pre- 
scribed degree distribution; 

2. Each node is labelled with a stylized banking balance sheet that identifies its 
external assets and liabilities, its internal (i.e. total interbank) assets and liabili- 
ties, and 7, its net worth or equity (i.e. its total assets minus its total liabilities). 
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Initially, the system is in equilibrium, meaning each node has positive net worth 
7 > 0. 

3. Each directed edge is labelled with a deterministic weight that represents the 
positive exposure of one bank to another. These weights depend deterministi- 
cally on the in-degree of the edge, and are consistent with the interbank (IB) 
assets and liabilities at each node; 

4. A random shock is applied to the balance sheets in the system that triggers the 
default or insolvency of a fixed fraction of nodes; 

5. The residual value available to creditors of a defaulted bank is zero, and thus 
the shock has the potential to trigger a cascade of further bank defaults. 

The principle of limited liability for banks means that equity holders are never asked to 
cover a negative net worth of an insolvent firm. Instead, the insolvent firm is assumed 
to "default" , meaning it ceases to operate as a going concern, and its creditors divide 
the residual value. Since this residual value is always less than the nominal liabilities, 
creditor banks thus receive a shock to their balance sheets, creating the potential for 
a default cascade. The GK model makes a very simple zero recovery assumption that 
residual values of defaulted banks will be zero, and thus every time a bank defaults 
a maximal possible shock will be transmitted to its creditors. 

The framework for random financial networks presented here generalizes the GK 
model in one important respect, namely that the edge degree distribution is arbitrary, 
allowing for any desired amount of assortativity or disassortativity in the network. 
Under these more general assumptions, we are able to prove our main result about 
default cascades, which draws the following conclusions: 

1. The probability that an edge is defaulted after n steps of the cascade depends 
on the in-degree j of the edge, but not on the out-degree. Let these quantities 
be denoted a^"^ = {a^"''*}jgz+5 ^+ := {0, 1,2,...} for each value of the in-degree 
j- 

2. The sequence {a^"'^}n=o,i,... satisfies a recursion 

a("+i) = G(a(")), n = 0,l,... 

where the monotonically increasing mapping G : [0, 1]^+ — >• [0, 1]^+ depends 
explicitly on the structure of the network and the initial shock distribution. 

An earlier result of this type was proved in Gleeson et al. [2011 in a situation in 
which the probabilities a are independent of j, and G reduces to a scalar function. 

One important implication of the above result is the so-called "cascade condition" 
which is a statement about the derivative of G computed in the limit as the initial 
shock goes to zero. It represents the condition that a tiny initial shock can create a 
finite size cascade, and is the main analytical quantity that determines the stability 
of the financial network. For example, with a given network structure, we are able to 
compute a critical value 7c of the net worth parameter that separates a no-cascade 
region from a cascade region. 
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The remainder of this paper is structured as follows. In Section [2] the extended 
GK model is described in detail, and the analytical description of its solution is 
given in Section [3j Section [4] discusses the cascade condition, and Section [5] derives a 
formula for the frequency of large scale cascades arising from an infinitesimally small 
seed. Numerical Monte Carlo simulation results and the corresponding analytical 
predictions are compared for a simple class of networks in Section [6j This section 
also investigates in more detail the relation between assortativity and systemic risk 
in a richer family of infinite network examples. Section [7] concludes. 

2 The Extended GK Model 

In this section we completely specify the extended GK modelling framework for which 
our analytical techniques will apply. The specification will consist of three levels: 
first, the random directed graph model for the "skeleton" of the network; second, a 
specification of balance sheet values corresponding to all nodes and edges; thirdly, 
a specification of the type of initial shocks that will be considered. We shall work 
within a probability space (JIat, J-at, IP) for some N < oo, where a general outcome 
a; G (7 AT is a directed graph with N nodes, with specified balance sheet values and an 
initial shock. 

2.1 The Skeleton Network 

The first step in building a random financial network is to build the skeleton directed 
graph that labels the banks and their interbank connections. Our construction is an 
extension of the Erdos-Renyi random graph model, and to describe it we introduce 
certain graph theoretic definitions and notation: 

1. A node v £ V = Uj^Vjk has type (j, k) means its in-degree is deg~(t;) = j and 
its out-degree is deg^{v) = k. We shall write k^ = k,jy = j for any v G Vjk and 
allow degrees to be any non-negative integer. 

2. An edge £ G £ = UkjCkj is said to have type {k,j) with in-degree j and out- 
degree k if it is an out-edge of a node with out-degree k and an in-edge of a 
node with in-degree j. We shall write deg~'"(^) = ki = k and deg~(^) = ji = j 
for any £ G Ckj ■ 

3. We write (or C~) for the set of out-edges (respectively, in-edges) of a given 
node V and (or vj) for the node for which £ is an out-edge (respectively, 
in-edge). 

4. We will always write etc. to refer to in-degrees while k, k' , k" , etc. refer 
to out-degrees. 

Figure [l] illustrates the neighbourhood of a type (3, 2) bank. Our random graph 
model is now characterized by a certain independence structure, and by probability 
laws P, Q for nodes and edges. 
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Figure 1: The skeleton structure of the network locahty of a bank v. Bank v is in the 
{j, k) = (3,2) class, since it has 3 debtors and 2 creditors in the interbank (IB) network. 

Definition 1. 1. For each j,k, Pjk := IP[Vjfc] is the probability of a type {j,k) 
node. This distribution has marginals := Pjk, Pj^ '■= "^Zk^jk- 

2. For each j, k, Q^j := P[£fcj] is the probability of a type {k,j) edge. This distri- 
bution has marginals := Yl,j Qkj,Qj ■= ^k Qkj- 

For any finite number of nodes A'^, a random graph with the required statistics can 
be constructed by following a generalization of the "configuration graph" construction 
used to construct Erdos-Renyi graphs: 

1. First, N degree pairs (j„, kn),n = 0,1, . . . , N are drawn independently from the 
P distribution. The nth node then has j„ "in-stubs" and k^ "out-stubs" . 

2. Next, for each in-stub attached to a node with degree pair {j,k),j > one 
draws an out-degree k' > from the conditional distributiorj^ F[ke = k\ji = 
j] = Qkj/Qj- we can say the type of this in-stub is {k',j). The collection of 
such draws is mutually independent. 

3. Then, each in-stub of type {k,j) in succession is glued to a random out-stub 
selected with uniform probability from all remaining unpaired out-stubs with 
out-degree k. Again, the collection of such draws is mutually independent. 

4. Finally, the graph is accepted if all in-stubs and out-stubs are paired, and re- 
jected if not. 

For a realization of the above finite graph construction to be accepted, it is 
necessary for the total number of in-stubs with any given out-degree k to equal the 
number of out-stubs connected to nodes with this value of /c. As — t- cxd, this leads by 

^Recall the definition of conditional probability for any pair of sets of outcomes A, B: 

V[A\B] ^V[AnB]/V[B] . 
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the law of large numbers to a consistency condition relating the marginals of Q to the 
marginals of P. Some further thought about the wiring of in-stubs to out-stubs leads 
one to deduce certain independence properties of the random graph in the N = oo 
limit. Heuristically, any node v or edge i divides the graph into "downstream" and 
"upstream" parts. We can say that any downstream probability is independent of 
the upstream graph and vice versa. Based on these considerations, we adopt the 
following description of the random directed graph model: 

Definition 2. For any N > I (including N = oo), let Pjk,Qkj,j,k = 0,1,... be 
collections of node and edge probabilities, subject to the constraints 

Q+ = kP+/z, Qj=jP7/z 

for all j,k, where the mean degree is z = "^^k^-^k ~ SjJ^j"- ^^''^ ^ ^ 
random directed graph model Q{N,P,Q) is defined by the above configuration graph 
construction. For N = oo, the random directed graph model Q{N = oo, P,Q) is 
instead characterized by the assumption of independence of downstream and upstream 
probabilities, conditioned on the type of node or edge. This independence condition is 
illustrated by the important cas^- For all j, j' , j" , k, k' , k" , 

P[v/ G Vjk,v^ G Vj'k'\i G Ck"j"] = Sj>jff6kk"^[jv = j\kv = k] ^ky = k'\jy = f] . (1) 

Remark 3. (Independent edge condition) The special case Q^j = kjP^Pj^/z"^ = 
QjQk '^^ss when each in-stub is glued to an out-stub selected uniformly from the 
collection of all un-paired out-stubs, and is the usual notion of random graph. We 
are interested in the more general case described above because observed financial 
networks do not appear to satisfy the independent edge condition. Apparently real fi- 
nancial networks have the "edge-disassortative property" that high degree banks attach 
preferentially to low degree banks. 

Remark 4. A natural measure of edge-assortativity by degree is the "edge-assortativity 
coefficient" rg G [—1,1] given by 

E,kJk[Qkj-QjQt] 

rQ = I ■ (2) 



This quantity has been measured for the Brazilian banking network in Cont et al. 



_ 2010^ and found to be usually negative. However, we will find some evidence that sys- 
temic risk of a network is related more strongly to a combination of edge-assortativity 
and node-assortativity (arising from the dependence between in- and out- degrees of 

A conditional independence structure more general than ([T]), not arising from the above graph 



struction, is analyzed in Boguna and Serrano 



2005 
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nodes). We therefore define a measure we call the " graph- assortativity coefficient" 
r G [— 1, 1] given by 



(3) 



where Bjj/ is the joint distribution of the in-degree of pairs of nodes connected by an 
edge: 

Bjji = P[j^ = j, jt,' = j'\v is joined by a single out-edge I to v'] 

BjkQkj' 



k 

and Bj = J2fBjj'^ B^, = J2j Bjj' are the marginals. Here, the formula for Bjj' is 
derived from Q . 

Remark 5. We observe that infinite graphs in our framework are essentially "acyclic", 
meaning the probability that any given node will be contained in a cycle of edges of 
any fixed length will be zero. This property is unlike what we observe in real financial 
networks. Since infinite graphs are almost everywhere tree-like, one can understand 
that the independence conditions ([T]) are part of a hierarchy of similar conditions for 
edges and nodes in general position. 



2.2 Balance Sheets 

To build a financial network with full accounting information, consistent with a given 
skeleton graph, one specifies the external assets Yy and external liabilities for each 
node V, and for each edge i of the network, an exposure size or weight W£. All these 
quantities are positive. From this one defines the net worth or buffer of a node v to 
be 

^y = Yy+ wi - Zy - Wi . (4) 

We will always assume that the system is initially in a "cascade equilibrium" (or 
"equilibrium" for short) in which all banks are solvent, which means that 7^ > at 
every node v. Thus can be thought of as a buffer that keeps the bank solvent when 
subjected to balance sheet shocks up to a certain size. 

The cascade dynamics in the GK framework do not depend on full accounting 
information, but only on the partial information 

{^y,v^v}l^{wl,^^C\. (5) 

We adopt a deterministic rule for which buffers may depend on the node type (j, /c) 
but the edge weights depend only on the in-degree deg~ {€) : 

( lv= 7jk, V G Vjk .g^ 

\ Wi = Wj, ie Ckj ' 
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For the original GK model described in Gai and Kapadia 2010 , which we will call 
the GK specification, the following choices are made: 



Ijk = 7 := 0-035; Wj = — , 

but the analytical results of that paper clearly hold for general prescriptions of the 
form ([g]). 



assets of bank / liabilities of bank / 







net worth yct^-y 






external assets 


deposits 








IB borrowings 
of bank / 




> 


IB loans by 
bank / 





Figure 2: Schematic balance sheet of banks in the (j, k) = (3, 2) class. 



2.3 Shocks and the Solvency Condition 

Insolvencies arise in a system initially in equilibrium only when a shock hits that is 
hard enough to cause at least one node to suffer a loss larger than its buffer 7^. For 
simplicity, we suppose that such an initial shock to our system causes an initial set 
Wo C V of nodes to become insolvent (for example by hitting their external assets), 
but leaves other banks' balance sheets unchanged. The set Wo is drawn randomly, 
with the fraction of type (j. A;) nodes that are defaulted denoted by 

P% ■■= nv G Wol^; G Vjk] . 

Under the GK "zero recovery" assumption that an insolvent bank can pay none of 
its interbank credit obligations, each insolvent node v G Wo triggers all its out-edges 
to have zero value. This triggering of edges to default is an instance of what we call 
an "edge update" step of the cascade: corresponding to any default node set W there 
is a default edge set C £ defined by the condition ^ G if and only if G W. 
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Each such defaulted edge I now transmits a maximal shock wi to the asset side 
of the balance sheet of its in- node (the creditor bank). If all balance sheets are 
determined by the reduced accounting information {7jfc} U {wj\^ then when is a 
set of defaulted edges, the solvency condition on a node v G Vjfc i^ 

This triggering of nodes to default we call a "node update" step of the cascade: 
corresponding to the default edge set M. there is a default node set W' defined by 
the condition v G W if and only if 

#{£- n Ai} < M,-fc := [TjfcMl , (7) 

where (j, A;) is the type of v. Here \x\ denotes the "ceiling" function, i.e. the smallest 
integer greater than or equal to x, and so Mj^ is the threshold for the number of 
defaulted in-edges that a type (j, k) node can sustain without itself defaulting. 

To summarize, for N < oo and N = oo, a specification {N,P,Q,'^,w, p^) of 
an extended GK financial system is the following information: (i) the "skeleton", a 
Q{N, P, Q) random directed graph defined by the probabilities Pjk, Qkj over all node 
and edge types; (ii) a reduced accounting set, denoted by {"fjk} U {wj} and (iii) an 
initial shocked set Wo with the default probabilities for each node type and the 
corresponding defaulted edge set A^o- 

Given any realization of an extended GK financial system so specified, the com- 
plete default cascade will be a deterministic alternating sequence of node and edge 
updates (finite if is finite), beginning with the initial shocked set Wo and its cor- 
responding edge set A^o- When = oo, the cascade can be fully resolved, and the 
expected fraction of total defaulted nodes and edges (and other statistics) can be 
determined by the inductive analysis given in the next section. 



3 Default Cascade Steps 

Given any realization of an extended GK financial system [N, P, Q, 7, w, p^) as spec- 
ified above, with an initial shocked set Wo and the corresponding edge set A^O; the 
default cascade can be thought of as a sequence of updates: 

/Wo\ /WoUWi\ /WoUW2\ /WoUW„\ 

VA^oy [moumJ \M0uM2J" [moumJ 

Note that the above set unions are assumed to be disjoint, since we prefer to distin- 
guish the initial default sets Wo,A^o from the sets of "newly defaulted" nodes and 

■^The indicator function 1^ of any set A is the random variable that is 1 on the set and on its 
complement. 
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edges. Inductively, we have increasing sequences of sets: 



yVn ■= defaulted nodes not in Wo "triggered" by edges in A^o U7W„_i (8) 
Mn '■= defaulted edges not in Mq "triggered" by nodes in Wo U W„. (9) 

When N = oo these random default sets Wn^-Mn define probabilities for n = 
0,1,2,... 

p% := F[vGWn\veVjk] (10) 

CT^,. := e Mn\i e Ckj] . (11) 

Now the set A^o is determined from Wo by an edge update step. Similarly, for 
each n > 1 the set Aio'J -Mn is determined from WoUWn by an edge update step. In 
all these cases the probabilities a'^j are determined by the following general lemma. 

Lemma 6. (Edge update) Let N = oo and suppose W C V denotes a set of defaulted 
nodes and for all j, k, pjk '■= IP[f € W|f G Vjk]- If the corresponding set of defaulted 
edges is denoted Ai C C then for all j, k 



Proof: When N = oo we may use the decomposition ([T]) and find 

j' 
j' 

= ^ Pj'k^ljv = j'\kv = k] 
j' 

Yuj'iPj'kPj'k) 



^k 



□ 



Given any set M. of defaulted edges, each of the nodes (excluding the initially 
defaulted nodes in Wo) recomputes its balance sheet and the node update step leads 
to a subset W C V \ Wo of defaulted nodes determined by M.. We separate out the 
originally defaulted nodes Wo since these were not triggered by defaulted edges. The 
probabilities associated to W are characterized by the following result. 

Lemma 7. (Node update) Let N = oo and suppose Ai C C denotes a set of defaulted 
edges with associated probabilities atj '■= IPi-^ G A^|£ G Ckj\- Then for all j, aj defined 
to be ¥[i G M\je = j] is given by 

tin 

Q- 
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If the corresponding subset ofV\ Wo of defaulted nodes triggered by Ai is denoted W 
then for all j, k, pjk := IP[f £ W\v € Vjk] is given by 

where the default thresholds are defined as in ^ by Mjf^ = \'yjf^/wj']. 
Proof: First we compute that 

„,=P[,eA<|i,=il = £4|^l^ = £i|p>. (12) 

Note that ¥[£ E M\v^ £ Vjk] = G M\je = j], that 

F[v e W\v G Vjfc] = F[v G G Vjk\m P[f ^ Wo] = {l-p%)nv e Wjr; G V.AWo] 

(13) 

and that a node v G Vjfc \ Wo wih be in default if and only if at least Mjk in-edges to 
V are in M. By the independence structure of the network for N = oo, the random 
variables Ii^m for all i G C~ , under the condition that v G Vjk \ Wo, are a collection 
of j identical independent Bernoulli random variables with probability aj. Putting 
these facts together gives 



which combined with (13) leads to the required result. 

□ 



Using these lemmas and the definitions (10), (11), it is straightforward to piece 
together the steps of the default cascade and obtain the main result of the paper. 

Proposition 8. Consider the infinite extended GK financial network (N = oo, P,Q,^,w, p^) . 
For n = 0,1, 2,..., let a'^"-'> = {a^"^} denote the probabilities F[£ gMqU Mn\je = j]- 
Then 

1. For n = we have 

ol = (14) 

(0) J2kiQkj(^kj) 
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2. For n = 1,2, . . . , the quantities p", o"", a^") satisfy the recursive formulas 

j 



Pjk 



m 



a 



kj 



(n) 



^k 



(16) 

(17) 
(18) 



where Mj^ = [Tjfc/t^jl- The total probability for defaulted {j,k) nodes at step 
n is p^jj. + and the total probability for defaulted {k,j) edges at step n is 

3. The new probabilities a^"^ = {oj"^} Q'^e a function G{a^"'~^^) which is explicit in 
terms of the specification {N, P, Q, c, w, p^) . 

Proof: In Part 1, (14) follows from Lemma [6] and (15) from Lemma [Tj Part 2 is 
immediate from the same two lemmas, while Part 3 is simply a composition of (16), 

(Tt]), 

□ 



4 The Cascade Condition 



The size of global cascades in an extended GK financial network with N = oo has 
essentially been reduced to solving the fixed point equation 

a = G{a) (19) 

by iteration of the mapping G. Scalar equations of this sort, giving the expected 
size of cascades on directed networks, have been previously derived in various con- 
Gleeson] t2008b] , |Amini et al.| (Mo) . In |Gleeson| |2008b| , the main focus is on 



texts 



while Amini et al. 





percolation-type phenomena (see also the undirected networks case |Gleeson|[2008al ), 
2010 considers more complicated dynamics but takes the limit 

where initial default fractions 



0. The case considered in Gleeson et al. 2011 



P 

can be difi'erent for each (j. A;) class, has not, to our knowledge, been considered pre- 
viously. In the current work, we include for the first time (through Qkj) the effect 
of non-trivial correlations between the degrees of nodes at either end of a randomly 
chosen edge. 



As a consequence of the Knaster-Tarski Theorem, equation (19) always has at 
least one solution and this will be a vector of probabilities a"^ S [0, 1] for all j. 
To see this one observes that G is a monotone mapping from the complete lattice 
[0, 1]^+ onto itself, under the partial ordering relation defined by a < 6 if and only if 
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cij < bj for all j G Z_|_, that is, G{a) < G{h) whenever a < b. This is enough to ensure 
the existence of at least one fixed point on the set [0, 1]^+. 

One important question is to consider initial points ea for small e > and ask 
whether the fixed points aoo(e) obtained this way are of order e or of order 1 as e — )• 0. 
In other words, what is the "cascade condition" that determines if an infinitesimally 
small seed fraction will grow to a large-scale cascade? This depends on the spectral 
radius of the derivative matrix D = {Djj/} with Djj' = 9Gj/5aj'|g=o,pO=o- Recalling 
that the spectral radius of a matrix D, \\D\\ := maxg.||^||^i 11-^^11; is the magnitude of 
the largest eigenvalue of D, one can see by a version of the Perron- Frobenius Theorem 
for non- negative infinite matrices, that the fixed points aoo(e) will be 0(e) if ||D|| < 1 
and will be 0(1) if > 1. Such a "cascade condition" plays a role in systemic risk 
analogous to the basic reproduction number Rq in epidemiology. In our framework, 
the derivatives Djji turn out to be easy to calculate and we find: 

Proposition 9. The infinite extended GK financial network (N = <x,P,Q,^,w, p^) 
satisfies the cascade condition, that is, any infinitesimal seed will trigger a large scale 
cascade almost surely, if the spectral radius \\D\\ > 1 where 

If \\D\\ < 1, then almost surely the network will not exhibit large scale cascades for 
any infinitesimal seed. 

In Section |6j we shall see that the cascade condition is indeed a strong measure 
of systemic risk in finite simulated networks. One can check that in the setting of 
independent edge probabilities, the mapping G reduces to a scalar function, a result 



that has been derived in a rather different fashion in Gai and Kapadia 2010 . They 



extend Watts' Watts 20021 percolation theory approach from his work on undirected 



networks to the case of directed nonassortative networks. We will see in the next 
section that the percolation approach to the cascade condition extends further to our 
directed assortative networks. 

We can understand the cascade condition more clearly by introducing the notion 
of vulnerable node, that is any node that becomes insolvent if any one of its debtors 
defaults. In our specifications, a (j, k) node is thus vulnerable if and only if < Wj. 
The matrix element Djj/ has a simple explanation that gives more intuition about the 
nature of the cascade condition: it is the expected number of edges £' with j> = j' 
that connect through a vulnerable node to an edge i with ji = j. Then for small 
values of a, one has a linear approximation for the change in a in a single cascade 
step: 

a^+^-af = Y,Djr {a])-a]l~') + 0{\\af) . (21) 
j' 

The condition for a global cascade starting from an infinitesimal seed is that the 
matrix D must have an expanding direction, i.e. an eigenvalue bigger than 1. 
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5 Frequency of global cascades and the giant 
vulnerable cluster 



The previous argument does not tell us directly about the frequency of global cas- 



cades. However, it is well-known [Newman 2010, Chapter 13.11] that the frequency 



of global cascades in infinite random graphs is given by the fractional size of the 
so-called in-component associated to the giant vulnerable cluster. 

To facilitate the discussion we make the following further definitions 

• V(^) is the set of vulnerable nodes. 

• Ss is the giant strongly connected set of vulnerable nodes (the "giant vulnerable 
cluster"); 

• Si is the set of (possibly not vulnerable) nodes that are forward connected to 
Ss by a path of vulnerable nodes (the "in-component" of the giant vulnerable 
cluster) ; 

• Tjk = l{7j.fe<«;,}- 

We are interested in the following probabilities b = (bk), bk ■= ^[v ^ Si\k^ = k] 
and note that v € Sf (i.e. the complement of Si) is equivalent to the condition that 
all the downstream nodes vj ,£ G are in the set {Sf (1 V(^,)) U V^'^^^. Thus, letting 
v' denote any node one edge downstream from v, one has: 

bk = (ck)'' , (22) 

where 

Ck = Yl e {Sf n V(,)) U Vl,)\v' G Vfk',ke = k]F[v' G Vj'k'\ke = k] . 
j',k' 



Next note 



[v' e {Sf n V(^)) U Vl^^)\v' G Vj>k', h = k\= Tyk'hk' + (1 - Tjvfc,) 

Pj'k'Qkj' 



[v G Vj'k'\ke = k] 



PfQt 



and thus 



Cfc = E i'^fk'bk' + (1 - r,,k')) ^J^^ . (23) 

Since b^ = {ck)^ it follows that c = (c^) satisfies the fixed point equation c = h{c) 
where for any sequence c = (c^) 

hk{c) = Y: (rMck'f + (1 - r,'.o) ■ (24) 



j',k' 



PfQt 
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Note that the equation c = h{c) has a trivial fixed point e = (1,1,...) that 
corresponds to the set Si having probabihty zero. We now verify that the cascade 
condition ||D|| > 1 is equivalent to the condition that e is an unstable fixed point, in 
which case there will be a nontrivial fixed point < Cqo < e . A sufficient (and almost 
necessary) condition for e to be an unstable fixed point is that ||D|| > 1 where the 
derivative D^k' = idhk/dck')\s=e is given by 

One can verify directly that 

D={ABAA-^)^, D = AB 

for matrices 

and from this it follows that the spectral radii and spectral norms of D and D are 
equal. Hence ||D|| > 1 if and only if ||D|| > 1. 

As long as the cascade condition is satisfied, the cascade frequency / equals the 
fractional size of the in-component Si and is given by 

f = Y^r[y(Z Si\k, = k] F[k, = k]= ^(1 - Ck,oo)P^ . (26) 

k k 

Repeating this type of argument to determine the size of the out-component of 
the giant vulnerable cluster, one also obtains an upper bound on the size of the global 
cascade. 



6 Numerical Results 

In this section we present results from large-scale Monte Carlo simulations on random 
networks, and show that the analytical theory of Section [3] for N = oo matches well 
to the numerical results when A^, the number of nodes in the network, is sufficiently 
large. 

6.1 A Simple Random Network Model 

We consider networks constructed with nodes of types (3, 3), (3, 12), (12, 3), (12, 12) 
and edges of the same types. We fix the marginal probabilities = P^ = 1/2 
which lead to an average degree z = 15/2 and the marginals = 1/5, — 4/5. 
For parameters a £ [0, 1/2] and b € [0, 1/5] the following P and Q probabilities are 
consistent: 

(P3,3 P3,i2A A/2-a a \ ( Q3,3 Q3,12\ ^ (l/5-b b \ 
\Pl2,3 Pl2,12j V « 1/2 - a) ' \Qi2,3 Ql2,12j \ b 4/5-bJ ' 

(27) 



15 



We first fix the value of a to be 0.5, which means that the in- and out-degrees of 
all nodes are negatively correlated: nodes with in-degree 3 have out-degree 12, and 
vice versa. We examine three different values of the parameter b: the independent 
connections case b = 0.16, the (almost) maximally assortative case b = 0.01 and the 
(almost) maximally disassortative case b = 0.19. Note that the independent edge 
condition has been assumed in all previous work on such problems. We also note 
that with b = 0, edges are maximally assortative and link nodes of out-degree 3 to 
nodes of in-degree 3 only, and nodes of out-degree 12 to nodes of in-degree 12 only. 
In this case, the network falls into two disjoint pieces. 

The balance sheet quantities are those of Gai and Kapadia |2010| (except for the 
percentage net worth 7, which we vary over the range 0% to 10%), while the initial 
shock distribution is taken to be p^j^ = 1/N for all types (j, k), corresponding to the 
shocking of a single, randomly-chosen, bank. 

Generating a finite size skeleton network with prescribed probabilities following 
the steps outlined in section 2.1 always encounters the difficulty of ensuring that the 
consistency condition for wiring the network is exact. The naive approach would be 
to accept any randomly drawn set of node types {{jn, kn)}n=i,...,N only if the consis- 
tency condition is satisfied: one can easily see that even for our simple probability 
specification the rejection rate will be extremely high, and the method very slow. 
Instead one in practice adopts some approximation that we call "clipping": given 
an inconsistent labelling, one selects random pairs of nodes, and "rewires" a small 
number of times to ensure the consistency conditions. Such a procedure may create 
a small bias that vanishes as — )• 00. 

Figure [s] compares theory curves for cascade size (found by iterating equations 
(14)-([T8| to convergence) as well as cascade frequency (given by (|26|)) with results 
from numerical simulations on random networks with N = 10^, 10 cind 10^ nodes. 
The nodal correlation parameter is fixed at a = 0.5, while the edge correlation pa- 
rameter takes the values b = 0.01,0.16,0.19. Results are plotted as functions of the 
percentage net worth parameter 7. In each case, 500 realizations are used to find the 



extent of global cascades (a global cascade is defined, similarly to Gai and Kapadia 



2010 , Gleeson et al. [2011| , as one in which more than 5% of nodes default), and 



the frequency with which such global cascades occur. As expected, the analytical 
approach accurately predicts the size of the global cascades. Some discrepancies may 
be noted in Figure |3j where the theory does not predict some global cascades, but 
note that these occur with only very small frequencies. 

The cascade condition ( 20 ) predicts that the critical values of the buffer parameter 
7 are: 7c = 0.067 for the parameters of Figure |3]^a), and 7c = 0.017 for the case of 
Figure . These values match very accurately to the locations of the dramatic 
transitions in the theory curve (and in the expected size of cascades in numerical 
experiments): for buffer values in excess of 7c global cascades are extremely rare, 
while for 7 values less than 7c the entire financial system is likely to fail following 
a single bank's default. These result indicate the potential usefulness of the cascade 
condition as a measure of systemic risk. 
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We consider in Figure|4]the joint dependences on a, b of various theoretical quanti- 
ties in the large limit. In the top figures, the critical value of 7 and cascade size are 
seen to be discontinuous, and not directly related to edge-assortativity (parametrized 
by b). On the other hand (see bottom figures), the frequency of cascades is contin- 
uously varying, and does appear to correlate somewhat with the graph assortativity 
coefficient r given by ([S]). 



6.2 Another Simple Random Network Model 



Now we have shown that the infinite N theory meets our expectations, we can further 
explore the implications of the analytical method. Since the specification of extended 
GK networks has many components, one must be rather careful in the questions one 
wishes to address: we choose here to try by means of a simple network specification 
to shed some additional light on the role the assortative properties of a network play 
in its susceptibility to systemic risk as measured by the cascade frequency. 

We consider stylized networks with many small banks and a few large banks. 
The set of node types will be {(2,2), (4,4), (8,8), (16,16)} with a diagonal node 
probability matrix: 

P := (Pjk) = diag(8, 4, 2, 1)/15 . 
The following edge probability matrices Q = {Qkj) are consistent with P: 
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and their convex combinations qiQi + q2Q2 + q.'iQ?,+ QaQa, 9i +92 + ^3 + Q'4 = 1, 
span a simplex of possible edge probability matrices. We can see that as measured by 
tq, Q(^) is maximally assortative, while Q^'^^ and Q^^^ are maximally disassortative, 
and the independent case is g(°) := [Q^^) + Q^^) + q(3) ^ q{4)]/4_ 

For the remaining components of the specification we adopt the default GK bal- 
ance sheet values but with 7 a variable parameter and consider shocking a single 
randomly selected node (this is an infinitesimal shock in infinite volume). We then 



compute the critical 7c using the cascade formula ( 20 ) , the cascade size from Propo- 
sition 8, and the default frequency using equation (26). 

Figure [5] shows how the theoretical values of 7c and cascade size depend on the 
particular Q matrix. Figure [6] shows how the theoretical values of the graph assor- 
tativity coefficient and cascade frequency depend on Q. In both figures, the four 
rows correspond to the simplices of Q matrices with 54 = 0, ^3 = 0, = 0, = 
respectively. 
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We see again in these networks that r and / vary continuously, while -jc and 
cascade size take on only discrete values. Of particular interest is the discernable 
covariation of / with r. Since r depends only on the skeleton graph and not on the 
balance sheet data, we cannot expect a one-to-one relationship between the quantities 
/ and r. However, we conjecture that r is in some sense the best possible purely 
graph theoretic measure of systemic susceptibility. Heuristically, we might expect 
that systemic risk is lowered if, all else being equal, the network is such that the 
correlation between in-degrees of neighbouring nodes is lowered. 



7 Concluding Remarks 

In summary, we have described here a rigorous analytical framework which can pre- 



dict the systemic risk of "deliberately simplified models" such as Gai and Kapadia 



2010 . The qualitative type of networks one can address has been extended con- 



siderably over existing work, in particular by the inclusion of the non-independent 



connections between nodes. The example of Subsection 6.1 demonstrates that finite 
size effects do not appear to dramatically impact systemic risk as long as > 100. 
More subtly, we also observed that graph-assortativity, rather than edge- or node- 
assortativity, can strongly affect the course of contagion cascades, and hence show 
the importance of incorporating assortativity in numerical and analytical treatments 
of banking network models. Our framework will enable extensive studies of alter- 
native network topologies; the cascade condition and cascade frequency provide two 
simple and useful measures of systemic risk by which to compare different network 
topologies. However, the daunting range of network variables means that both ana- 
lytical and numerical studies must be carefully framed to address specific issues, for 
example, to uncover other key determinants of systemic risk. Finally, we anticipate 
that future work can show how the approach described here may be further extended 
to include partial recovery models (such as Nier et al. [2007] ) and stochastic balance 
sheets. 
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Figure 3: Numerical simulation results (symbols) and theoretical results (curves) for the 
random network model of Equation (27), on networks of nodes with parameter a = 0.5, 
as functions of the net worth 7. The average size and frequency of global cascades in 
simulations are shown by red circles and blue crosses, respectively. Theoretical results for 
the expected cascade size (black solid curve) are from Section 3; those for the frequency 
of cascades (dashed magenta curve) are from Section 5. Each column shows results for a 
different network size A^, and the parameter b takes a different value on each row of the 
figure. 
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Figure 4: Theoretical joint dependences on the graph parameters (a, b). Top left: critical 7 



value for the random network model of Equation (27). Top right: Expected size of cascades 
(from Section 3) when 7 = 0.05 and p'jj^ = 10~^. Bottom left: the graph assortativity 
parameter r. Bottom right: Frequency of cascades (from Section 5) when 7 = 0.05. 
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Figure 5: Critical value of 7 (left column) and expected size (right column) for cascades on 
the random network model of Section 6.2 with 7 = 0.0375 and p^j^, = 10~^. The triangles 
shown correspond to Q matrices with (from top to bottom) = 0, gs = 0, q2 = 0, and 

gi = 0. 
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Figure 6: Graph assortativity parameter r (left column) and frequency (right column) for 
the same parameters as Fig. [5j 
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